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Abstract We address the question of representativeness of a single long unstable 
periodic orbit for properties of the chaotic attractor it is embedded in. Y. Saiki and 
M. Yamada [Phys.Rev.E 79, 015201(R) (2009)] have recently suggested the hy- 
pothesis that there exist a limit distribution of averages over unstable periodic orbits 
with given number of loops, N, which is not a Dirac <5-function for infinitely long 
orbits. In this paper we show that the limit distribution is actually a 5-function and 
standard deviations decay as N~ x l 2 for large enoughs. 



1 Introduction 

Recent investigations lfTTi r7ll6l fT2l[T0l have arose the question of representativeness 
of a single unstable periodic orbit (UPO) for properties of the chaotic attractor it 
is embedded in. Interestingly, while Refs. Q [6] [T2) concern unstable time-periodic 
orbits in turbulence, Ref. [10] discusses turbulent pipe flows in relation to spatially 
periodic solutions. Motivated by Q [321 [TOl Y. Saiki and M. Yamada flj] discus 
the calculation of average values along UPOs (unstable periodic orbits embedded 
into the chaotic attractor) and using the results of such a calculation for estima- 
tion of averages along a chaotic trajectory. In particular, with simple paradigmatic 
chaotic systems — the Lorenz system, the Rossler one, and a 6-dimensional eco- 
nomic model — Y. Saiki and M. Yamada explored the convergence of average values 
calculated along UPOs as the number of loops, N, of a UPO grows. They surpris- 
ingly found that the distribution of average values along all the UPOs with given N 
(henceforth, A^-UPOs) practically does not shrink for long orbits (N = 10. ..14 for 
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Fig. 1 Chaos arising via bifurcations of homoclinic orbits. The phase portrait and the Poincare 
section for the Lorenz system with classical parameter values (a). The Poincare map for (x+y) (b) 
is qualitatively similar to the saw map (c). 



the Lorenz system, N = 10. ..17 for the Rossler one), suggesting the existence of 
a limit distribution of averages over A^-UPOs, which is not a Dirac 5-function for 
infinitely long orbits (N — > °°). In Ref. lfl3l . it was shown that the UPOs with the 
length N <14 are not enough long for conclusions on asymptotic behavior. 

We are interested in the asymptotic behavior of the distribution of average values 
over UPOs with fixed N. Notice, the evaluation of average values for the chaotic 
attractor requires this distribution to be weighted with the distribution of the natural 
measure over the set of UPOs; in Ref. flU the natural measure of an UPO was shown 
to be inverse proportional to its multiplier. Meanwhile, we restrict our consideration 
to the question of convergence of the limit distribution and do not consider the dis- 
tribution of the natural measure. 

In this work we treat two maps: "saw" map 

x n +i=2x„ modi (1) 

— paradigmatic model for the chaos (i) arising via cascade of bifurcations of homo- 
clinic orbits (as in the Lorenz system JU, Fig.QJ) 10 19) — and "tent" map 

!2x n for x n < 1/2, 
(2) 
2 — 2x n for x„ > 1/2 
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Fig. 2 Chaos arising via cascade of the period-doubling bifurcations. The Poincare section (a) and 
map (b) for the Rossler system with classical parameter values . 

— paradigmatic model for the chaos (ii) arising via cascade of the period-doubling 
bifurcations (as in the Rossler system, Fig.|2^) 0131 [TJ. The reason for our choice is 
not only paradigmaticity of these maps but also the fact that one gains exceptional 
technical opportunities for calculating all UPOs in them for giant N. We find the 
distributions of averages upon A-UPOs to shrink as A grows; standard deviations 
decay as 1 / y/N. Then we consider the Lorenz (the first sort of chaotic systems) 
and find for it the same kind dependence of standard deviations on A as for the 
paradigmatic map models. Convergence of the distribution to a 5-function can be 
observed for enough long UPOs (for the Lorenz system A > 10, cf. Fig.HJi. 



2 Limit Distribution of Averages 

Let us first recall the origin of maps (Q~|i and (0 in order to make the physical repre- 
sentativeness of these models evident. In the systems where the both separatrices of 
the saddle knot P come back in the vicinity of the stable manifold W s of P (Fig.QJ) 
the touching of a separatrix and W s results in appearance of the homoclinic orbit. 
The cascade of bifurcations of such orbits leads to the formation of the chaotic set 
(in the Lorenz system this cascade occurs at the single point r « 13.927 owing to 
the symmetry (x,y) O (— x, —y), cf. J8][5l[9)). For the Lorenz system 




dy 

— = rx — y — xz, 




(3) 



with the classical parameter set [a = 10, b = 8/3, r = 28] we choose z = [(r — 
l)/(4b)] l / 2 \x + y\ as a Poincare section, and find that the transversal structure of its 
intersection with the chaotic set is very narrow (Fig.Q]i; the Poincere map for (x+y) 
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Fig. 3 The number of UPOs used for averaging is presented by circles for the saw (a) and the 
tent (b) map; until the plateau, all UPOs are detected and employed, and on the plateau we utilize 
UPOs chosen randomly with a probability uniform over all Af-UPOs. The standard deviations of 
{(x— 1/2) 2 ) for the saw map (c) and the one of (x) for the tent map (d) decay as \/\/N for enough 
large N. The PDFs of the averages can be seen to shrink slowly as N grows for both the saw (e) 
and the tent (f) maps; circles: N = 1 1, squares: N = 23, diamonds: N = 61. 



is visually unambiguous (Fig.QJi) and similar to the saw map (Fig.Q]:). The averag- 
ing of a certain value along a certain trajectory on the chaotic set can be reduced to 
the averaging of a certain function /(x+y) along the trajectory in the Poincare map 
of (x + y). The system symmetry suggests an even function f(x + y) = f(—x — y). 
Actually, the contribution of /(x+y) into (/(x + y)) over an orbit after one Poincare 
recurrence is weighted by T (x + y) / (T (x + y)) map , the ratio of the recurrence time 
for the trajectory running from (x + y) and the average recurrence time for the or- 
bit, (...) ma p stands for averaging over iterations of the Poincare map (not over real 
time). Hence, (/(x + y)) = (f(x + y)T(x+y)) ma? /(T(x + y)) ma? , i.e., (/(x+y)) is 
subject to the additional dispersion due to nonuniformity of (T(x + y)) ma p over var- 
ious orbits. However, the nature of the dispersion of (r)map 

is the same as the one 
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Fig. 4 (a): STD normalized by STD for N = 15 is plotted by circles for the saw map and by 
crosses for the Lorenz system [STD((x 2 )) 15 _ UP o s ~ 0.5086323] (b): The PDF of (x 2 ) for the 
Lorenz system shrinks only slightly as N grows at the interval [10,20]. Solid line: N = 11, dashed 
line: N = 15, dash-dotted line: N = 19. 

of (/r) map , and in this paper we do not introduce any additional complication into 
our averagings with the paradigmatic map models. On the whole, the possibility to 
construct the Poincare map in the vicinity of the saddle knot qualitatively similar to 
the saw map is the common peculiarity of the systems with a stable chaotic set of 
the kind we consider. 

With regard to the second sort of chaotic systems, the existence of a unimodal 
map is a key feature needed for the cascade of the period-doubling bifurcations to 
occur resulting in chaotic behavior of the system |2][3][T]. Thus, for instance, the 
Poincare map of the Rossler system (Fig. [2) is quite similar to parabola of the logis- 
tic map x„ + i = bx n (l —x n ). The attracting chaotic set of the logistic map possesses 
its largest size at b = 4 where this map can be turned into the tent one (f2]i by virtue of 
substitution x — > sin nx. Hence, the tent map is quite representative for the systems 
where chaos arises via the cascade of the period-doubling bifurcations. 

As it is noticed above, the symmetry of the Lorenz system suggests us consider- 
ation of (f(x)}: f(l /2 +x) = f(l/2 — x) along trajectories in the saw map. Hence, 
for simplicity, we consider ((x— 1/2) 2 ) for the saw map (Q]) and (x) for the tent 
one (Q~|i. For the Lorenz system we calculate (xr) and (z). 

The dynamics of the saw map ([]} can be easily dealt with within the frame- 
works of the binary notation of x. In this notation, an iteration of map ([TJ results 
in the shift of the binary point in x by one position to the right and omitting the 
integer part. Thus, x with an -periodic sequence of digits in the binary mantissa 
belongs to an Af -periodic orbit. Employing this fact one can strictly calculate giant 
amount of UPOs and averages along them. Dealing in such a fashion with the tent 
map d2j is less efficient and more sophisticated, but still possible. In Fig.|3]one can 
see that for enough large A the standard deviations of the averages decay as 1 / \/N. 
Additionally, for (x) in the saw map, which is actually beyond our immediate in- 
terest, one can analytically find probability P((x) = m/N) ps Nl/[2"ml(N — m) !], 
where m = 1,2,..., A— 1, and, for A > 1, P((x)) (27rA)" 1 / 2 exp(-2A(x) 2 ) and 
STD((*)V_ UP os = (4A)- 1 / 2 + 0(A- 1 ). 

In order to verify how far our findings for the paradigmatic map models are rel- 
evant for the chaotic systems with continuous time, UPOs in the Lorenz system (O 
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have been numerically calculated with double precision. We have calculated all 
UPOs with N < 19 (in Ref. JTT| 90% of UPOs with N < 14 were detected) and 
evaluated (x 2 ) and (z) along the detected UPOs. In Fig.|4]one can see the same scal- 
ing law for long UPOs; STD N~ x l 2 . The dependence of STD on N for the Lorenz 
system and the saw map are similar for small N as well (Fig.Hk). 

Notice, the tendency of the PDF to shrinking as N grows can be hardly detected 
from plots of PDFs (see Figs. |3^,f and|4j5) for N £ [10,20] which were typical for 
Ref. 03]. 



3 Conclusion 

Summarizing, we have analyzed two paradigmatic map models: the saw map ([TJ 
representing chaos arising via cascade of bifurcations of homoclinic orbits and the 
tent one (O representing chaos arising via cascade of the period-doubling bifurca- 
tions. For the both models we have reliably established the fact that the distributions 
of averages along unstable periodic orbits with given number of loops N shrinks to 
a Dirac 5-function for N — s- °°; the standard deviations obey the decay law N~ 1//2 for 
N >\Q. For the Lorenz system the same features have been confirmed. In particular, 
the average over a long UPO gives a correct representation of the average over the 
whole chaotic set which it is embedded in. 
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